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Abstract 

Using the Bjorken expansion model we study the effect of intrinsic hydrodynamic noise on the 
correlations observed in heavy-ion collisions. 

Rich data on fluctuations and correlations in heavy-ion collisions raise many 

interesting questions, the answers to which could shed light on the properties of the medium 
created in these collisions. For example, it has been understood recently that long-range rapidity- 
independent correlations are caused by the fluctuations of the initial collision geometry which 
are propagated hydrodynamically |5]. We wish to highlight the fact that in addition to the initial 
state fluctuations, which are outside of the scope of hydrodynamics, there is an intrinsic source of 
fluctuations in hydrodynamics itself — local thermal noise. We explore the eff'ect of this noise on 
observed correlations in a longitudinally expanding fireball created in heavy-ion collisions |6]. 

What is the origin of the intrinsic hydrodynamic noise? The local thermal equilibrium that 
underlies the hydrodynamic description is a statistical concept. The thermodynamic state is not a 
single state, but an ensemble of states, macroscopically similar, but with microscopic differences 
fluctuating from one member of ensemble to another member For example, the equation of state 
expressing pressure as a function of energy density, P - Pie), gives the average value of the 
pressure. The pressure fluctuates even if the energy density is fixed (picture a gas of molecules in 
a box hitting its wall). Had there been no fluctuations, correlators such as, e.g., ( T/jyT^.^ ) (which 
give viscosity via Kubo formula), would have vanished. 

The relationship between fluctuations and dissipation suggests the idea that measurements 
of fluctuations may yield information about dissipative coefficients, such as viscosity, comple- 
mentary to existing methods based on the azimuthal asymmetry of flow. 

To see how the noise enters in hydrodynamic equations it is helpful to consider the origins of 
hydrodynamics. We know that a perturbed fluid relaxes back to an equilibrium. There are two 
very well (parametrically) separated scales associated with this process. First, the local thermal 
equilibrium is established on a relatively short time scale characteristic of microscopic processes, 
such as collisions. Achieving the global equilibrium, i.e., the same thermodynamic conditions 
throughout the volume of the fluid, takes considerably longer time, which grows as the square of 
the typical size of non-homogeneities. Hydrodynamics describes this much slower process. 

Hydrodynamics is an effective theory which deals only with degrees of freedom that matter - 
macroscopic quantities characterizing the local thermal equilibrium such as energy, momentum 
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and charge densities, which change slowly because the corresponding quantities are conserved. 
Faster microscopic degrees of freedom left out of the hydrodynamic description are the noise. 

To formalize the above picture we follow the approach of Ivf] and apply it to a relativistically 
covariant formulation of hydrodynamics. In this formulation, the covariant degrees of freedom 
are the components of the 4-velocity of the local rest frame (the frame where the momentum 
density r"' vanishes) and the energy density s in this frame. We can express 4 independent 
components of T'^'' in terms of the variables s and u^'. 

To close the system of equations the remaining 6 components of the T'^^ (the stress in the 
local rest frame) must be also expressed in terms of s and m''. In equilibrium T^'^ is constant 
throughout the system and Lorentz invariance together with definitions of s and m'^ mean that the 
stress tensor must be given by 

T^^ ^ eu^'u'' + P(e)A''\ (1) 

where Ay = 6y + u^^Uy is the projector on the spatial hyperplane in the local rest frame (A^m'' = 0). 
The deviations from the equilibrium are due to gradients and to the first order there are two 
covariant expressions, corresponding to shear and bulk viscosities: 



■ ^A''''(V ■ u). 



(2) 



The main point is that the expression T'"' = T'Z + Ar'"' holds only on average. Both sides of 



eq 

this relation fluctuate. For every member of the ensemble there is a discrepancy 5^^": 

vise 



(3) 



The discrepancy is the noise from the fast microscopic modes and, therefore, the correlation 
functions of this noise are local on the hydrodynamic scale: 



{s''\x)S"\y))-5\x-y). 



(4) 



The magnitude is determined by the condition that the equilibrium probability distribution 
of B is given by the number of microscopic states, i.e., exponential of the entropy at that e. In 
particular, dissipation due to shear and bulk viscosities must be matched by noise 



( S^'\x)S "^(y) )^2T J] [a^^A^P + A^/^A"") + - ^^j A'^'A''^ 



(5) 



Although the noise is local, the correlations induced by it are propagated by hydrodynamic 
modes over macroscopic distances. Stochastic equations 
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(6) 



solved for fluctuations of e and m'' around a static equilibrium state give well-known equilibrium 
correlation functions (e.g., Kubo relation for viscosity). Our goal is to apply these equations to 
determine correlations in a non-equilibrium, expanding fireball. 

The simplicity and symmetry of the Bjorken solution allows analytical treatment. Further- 
more, in this work we only consider rapidity dependence (i.e., we integrate/average fluctuations 
over transverse directions). 

We write the hydrodynamic equations in Bjorken coordinates and linearize in the fluctuations 
of e and u^. The typical structure of a correlation function at Bjorken freezeout time, Tf, is given 
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by the integral over Bjorken rapidity ^ and the Bjorken time r reflecting the fact that the noise, 
which is the source of the correlations, exists at all space-time points: 



<p(fl,Tf)p(^2,Tf)) = - — ^ fl'^Gp(^l-^;Tf,T)Gpfe-^;Tf,T) 



(7) 



where we use convenient notation p = 6s I s for the relative entropy density fluctuation and 
V = (4/7/3 + 01 s for "longitudinal" viscosity to entropy ratio. The fluctuation p at point 
is sourced by the noise fluctuation at each point ^ at earlier time r propagated hydrodynamically 
via the linear response Green's function Gp(^i - ^; Tf , r). The same fluctuation is also propagated 
to point ^2 via Gp(^2 - ^; Tf , t) and the correlation is the averaged product of the two fluctuations. 

In order to demonstrate explicitly the contribution of each Bjorken time slice r to the final 
correlation function we consider separately the expression in square brackets in Eq. (|7ll, which we 
denote as Gpp(A^; Tf, r). If one neglects viscosity, the correlator Gpp{A^; Tf, r) will have singu- 
larities at the separation equal to twice the distance traveled by sound on top of the expanding 
medium, namely, 2vsln(rf/T). This singular contribution is easy to understand since it results 
from a noise fluctuation located exactly halfway between the points and ^2 where correla- 
tion is measured. We separate the singular contribution (defined as linear combination of step 
function and its derivatives) and the remaining regular (more precisely, continuous) contribution 



Gpp = Gpp^ + G'pp and plot them using characteristic values of parameters in Fig.\T\left,center). 






Figure 1: Left: Singular pait of the contribution of a single time slice t. Vj^ = 1/3 and ln(Tf/T) = 4. In order to 
visualize the infinitesimally naiTow singularities we "smear" the function with a Gaussian of square width cr^ = 0.1. 
Center: The regular part of the correlation from 3 time slices t: ln(Tf/r) = 2,4,6. Right: Rapidity gap dependence of 
the particle-particle fluctuations in Eq. (9)- 



If the dispersion (jj{k) of small fluctuations around the Bjorken solution were linear, there 
would only be the sound-front singularities at - 2vsln(Tf/T) and 0. However, u - ia + 
sjvs^k^ - cP- \a = (1 - Vs^)/2]. Note that the slowest mode behaves as w_ ~ -ik^ for ^ 0. 
Correspondingly, the "wake" in Fig. ^center) approaches a Gaussian for large time intervals 
ln(T//r). Although diffusion-like in appearance, this mode is present even without dissipation. 

The main effect of viscosity is, as usual, to smear the sound-front singularity with a Gaussian 
of width which grows with time interval between Tf and t (initially following diffusive law cr ~ 
y/Tf - r): 



1 

TfTf 



(8) 



Finally, we translate fluctuations of hydrodynamic variables to observable fluctuations of 
particle distributions dN/drj using the Cooper- Frye freezeout formahsm. It takes into account the 
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fact that particles within a given local subvolume at Bjorken rapidity ^ have kinematic rapidities 
rj which are thermally spread around ^. The result can be presented in the following form: 



dN dN\ldN\ 45c/s 
6—d- \' \ 



di]i dT]2 l\ dr] I Ajr'^NefiiTo) TfT{ 



7-2 ^I's "-2 
^ 



KiM), (9) 



where we separated the dependence on rapidity gap At] - rji - 772 into a factor K(Arj). This factor 
does not depend on other parameters such as initial time and temperature tq, Tq, the effective 
number of degrees of freedom defined via s{T) - 27r^Neff{T)T^ /45 and the mass degeneracy of 
the observed species t/s (e.g., 2 for charged pions). It does depend on the ratio of the particle 
mass and the freezeout temperature and for pions at Tf = 150 MeV it is plotted in Fig. right). 

We conclude by describing a qualitative picture of the two-particle correlations as a function 
of rapidity separation which may be helpful in disentangling the two major contributions to 
correlations: initial state fluctuations and hydrodynamic noise. The initial state fluctuations, 
being fluctuations of the two-dimensional density of the colliding Lorentz-contracted nuclear 
"pancakes", create correlations which are largely independent of Arj and stretch over the whole 
rapidity interval between spectator fragments. In contrast, hydrodynamic fluctuations are local, 
and it takes time for them to spread over larger A77. The effect of the hydrodynamic noise is 
strongest at small Arj < I leading to a characteristic thermal peak (see Fig. right)). At larger 
At] the hydrodynamic fluctuations induce correlations whose strength decreases with Arj and 
vanishes beyond the sound horizon at 2vs ln(Tf/To). At such large A77 the correlation should 
reach a plateau set by initial state fluctuations. 

The magnitude of the noise-induced long-range correlations is proportional to viscosity, 
Eq. (|9]l. This fact could potentially be used to measure or constrain the value of this important 
transport coefficient by studying the A77 dependence of the correlations. 

One can extend this analysis to the dependence of the correlations on both rapidity and az- 
imuthal angle separation (Arj and A(p) in a way similar to the A0-only dependence analysis in 
Ref.|0]. This work is in progress (see T. Springer's contribution to these proceedings, |9]). Ex- 
perimental data on Arj, A(f) dependence of fluctuations, or, upon Fourier transform with respect to 
A(p, the Arj dependence of v„„(Ari) could be compared to such a theoretical analysis. Of course, 
a more realistic (e.g., 3d event-by-event hydro |10]) simulation may be necessary to enable a 
quantitative comparison and to determine or constrain the value of the viscosity. 
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